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SIMtARY 


Flutter teats were made of two wind-tunnel-fan aiodels, one of 
which had conventioiial Clai-k Y airfoil soctions end the other of 
which had high-coiiiber blade secticno. The restilte of the tests 
confirm the stall -flutter theory which predicts much higher flutter 
Bijeeds at the high lift coefiicionta for the hlgl'i-camber blades. 

Tiae higli-camber blades could, therefoi’e, be operated at much 
higlier power leading than could the low-camber blades. Aerodynamic 
data of these tests indicated little difference in efficiency for 
the two fans but somewhat highez' maximum lift coefjTlcient for the 
high-ceuzber fan. The efficiency of the fans decreased ai’ter the 
velocity, as calculated by two-dimensional theory, exceeded tJio speed 
of sound on the upper blade siu’faca . 

An analysis is JTiade of the factors that determine the power 
loading of a fan blade section. Graphs are presented which give 
the maximum power loading foi’ idealized sections of various 
thiclaiess I’atios operating at the critical Mach number. 

Examples show that the ideal oi- design lift coefficient of 
an airfoil is s-lmost the same as the lift coefficient giving the 
maximum flutter speed . It is tlierefore desirable that a section 
be operated at the ideal lift in order to obtain high critical speeds 
as well as the maximum margin of safety with respect to flutter. 


INTRODUCTION 


An Important problem in the design of ducted fans, such as 
wind-tunnel fans end axial comproosors, is the absorption of maximum 
power at high efficiency without blade failure. Tlie present paper 
is concerned with two of the factors, which may limit the po\ror 
loading of fans. These factors are the flutter speed and the 
critical Mach nuiibor of the blade section. 



2 


IIACA TIT No. 1330 


The flutter speed is shown in reference 1 to depend on the 
lift coefficient of the blade. For operation at hi^i lift coeffi- 
cients, theory shows that the flutter speed can be increased by 
properlj'' cambering the blade section. In order to check this 
theory, a model fein having high-camber blades was built and 
tested and the results were compared with those for a fan having 
Clark Y airfoil sections. The results of these tests are reported 
in the present paper. 

In order to determine to what extent hl^ lift coefficients 
can be used to advantage, an analysis is made in which the critical 
Mach numbers of the blade sections are teiken into consideration. 
Relations for the power loading and the pressure rise are obtained, 
and the maximum power loading per unit blade area is given for 
idealized airfoil sections of different thickness ratios operating at 
the critical Mach number. The analysis applies only to a blade 
element operating at ideal conditions and, therefore, cannot 
be directly applied to the perforBiance of the entire blade. The 
analysis serves, however, to show the important parameters and gives 
an upper limit to the useful power that can be absorbed if the 
section critical Mach number is considered to correspond to the 
limiting speed for efficient operation. 


SIMBOLS 


T thrust, pounds 

V axial velocity of stream, feet per secoM 

p density, slugs per cubic foot 

V resultant velocity at blade section, feet per second 

c chord, feet 

A blade area, sq.uare feet 

M Mach number 

M^ helical tip Mach number 

a speed of sound in air, feet per second 

p power, foot pounds per second 

power loading per unit blade area, horsepower per sq.uare foot 
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lift coefficient 

<if geometric helix angle, degrees 

r radius to propeller section, feet 

E propeller tip radius, feet 

B number of blades 

o solidity 

Cp presswe coefficient 

slope of lift curve 

^diq/ 4 pitching-moment coefficient about q[uarter -chord point 
X distance along x-axis, chords 

y distance along y-axis, chorda 

^c.g. location of center of gravity as measured from leading 
edge, chords 

Cr lift coefficient for ideal no -twist condition 
I 


untwisted or design value of lift coefficient 

5 dynamic pressure of operating speedy poxsads per square 

foot 

K torsional stiffness of blade, foot-pound per radian 

G shear modulus of elasticity, pounds per square foot 

J torsion modulus of section, feet^ 

t thicloiess of section, feet 

L representative length of blade, feet 

Q quantity rate of flow throu^ fan, cubic feet per second 

Ap pressure rise through fan, pounds per square foot 


^cal calculated classical -flutter tip Mach number 
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Mf Biaxiinvan measured flutter tip I-iach nuiater 

•^max 

Mf> mlnimuTd neaoured stall -flutter tip fiach nuidter 

rain 

a section angle of attack, decpreea 

P tlade angle at radius r, degL'eea 

angular velocity of propeller, radians per second 
p efficiency 

a' rotational -velocity interference factor 

Sutacripts : 

u untwisted or desltin 

cr critical 

I ideal 

div divergence speed 

f flutter speed 

i incompressille 

M conditions at critical Ma,ch niffiiter 

max raaxi-uum 

0 standard conditions 

THEORETICAL MALYOIS 

Relations for 2ovrer Tjcading and Pressure Rise 

The following analysis is hased on simple hlade -element theor^^ 
and neglects the drag forces and the change of densitjr of the fluid 
passing thi’ougli the fan. The quantities considered are shown in 
figure 1 . Tlie usefxil power is given by tiao T)i*oduct of the thrust T 
and the axial velocity V. 
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For oji elerient of "blade radius dr the relatione fc;r the 
power are as follows: 


dl^ = V dT = V dl. cos (f> 


where 

dL = i-pl/^O-rC dr 
V = \l sin 

Therefore , 


d? = cos 


^ cin 5^ c dr 


Since 


c dr = d/. 


end 


cos sin sin 2'^ 

“i" 


( 1 ) 


The quantity dp/dA is the power loadint^ per unit "blad.e ai'ea. 
When the resultant "blade -eleiiient velocity W is expressed in 
terns of Mach nini"ber M and speed of soimd a, equation (l) nay 
"be written 


df' 




sin 2^ 


If the geometric helix angle is 45^, the power loading per unit 
"blade area is a maxirnumj that is, 
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A convenient ©xpre salon for the maximum pover loading per unit 
"blade area Pa expressed in hcrsepovTer per siiuare foot for air 

at standard conditions is 


"0.002378 X 11163 


^^max 


550 


ClM- 


= 1500 


(h) 


Equation (4) may "be writtoin in a moi’s general fom for any density 
and speed of sound as 


P 


Amax 



Ci;i3 


( 5 ) 


where Pq refers to density of air at stand8j:‘d conditions aM 
a^, to speed of sound in air at standard conditions. 

An expression for the pressure rise may "be o"btained "by equating 
the tlirust of the fan "b-lade elements to tlie product of the pressure 
rise Ap times the area through which the "blade elements s^reep. 
Thus, for a fan having a given number of blades B and. a blade - 
element radius dr 


B dT = Ap2itr dr 


Ap = 


B dT 
2nr dr 


~pV?^CxjBc dr cos ^ 


2itr dr 
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Since 


Be 

2jtr 


a 


Ap s cos 5^ 


C6) 


Or, In terms of M and a, 


Ap = -pCjM^a^c C09 ^ 


(t) 


It may be noted that the maximum value of Ap occurs when ^ 
is 0°; whereas the maximum value of occurs when 0 is 45°. 


Maximum Power Loading as Limited 
by Critical Mach Number 

Equation ( 3 ) shows that the power loading is proportional 
to CjM°. Increasing the lift coefficient generally decreases the 

critical Mach number. The problem Is therefore to find the value 
of Cj^ and the associated section ci’ltical Mach number M^j, such 

that the product of is a maximum. In the present analysis 

the critical Mach number of the blade section is assumed to be the 
upper limit of efficient operation. 


Idealized secti ons.— Max. A. Heaslet of the Ames Ae;ronautlcal 

Laboratory of the NACA has sboTm the relation between critical Mach 
number and lift coefficient for various airfoils. As a limiting 
case an idealized section with elliptical thickness distribution 
was used; this section carried lift wl.th a flat— top lift distribution 
in which one— half of the lift was considered to act on the lower 
surface and one-half, on the upper surface. Such an idealized section 
has been used in the present paper to evaluate the quantity 
for various thickness ratios. 
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The method usod for calculating the critical speed of an 
airfoil at a given lift coefficient consist& of first calculating 
the potential tvo -dimensional pressure' distrihution. (See references 2 
to A.) The peak negative pressure is then corrected for compressi- 
hillty, and the stream Mach numher is found at which the local Mach 
nualber is unity. The lift coefficient for incoiipressihle flow is 
then corrected for compressihility to this stream, Mach number. This 
lift coefficient is designated Ci^j and the corre .spending stream 

Mach number is designated M^j.. For simplicity and consistency 

"both the peak pressure coefficient and the lift coefficient have 
"been corrected in the present paper hy Dieans of the Prandtl-Glauert 
f actoi’ , as follows : 


and 






/l - Mcr^ 




ll. 


M . /. 


j/l - Mcr*^ 


3 

In figure 2 the values of Cp^^cr "^•xij'num pov/er loading 

P. (hp/aq, ft) are plotted against for tlie fairiily of 

Aiiax 

hypothetical idealized airfoil sections- Lines of constant critical 
Mach number are also shomi. The figure shows tliat for these 
Idealized airfoil sections, the power loading becomes greater for 
the higher lift coefficients. For exaaple, a 9~percent~thick airfoil 
section operating at =0-3 has a critical Mach number of 0.8 

and a power loading of 225 horsepower per sq_uare foot. A section 
of the sajiie thiclcness operating at = 1.05 has a critical 

Mach nvimber of 0.7 and a power lo.ed.ing of 5^ horsepower per 
square foot. The maximum power loading occurs at the higliest lift 
coefficient at which the section operates. If the curves were 
extended far enou{yi, they would have maxisium values of 

the curves were calculated only to a value of 1.6 for sj.nce 

this value already is hij^ier than that norrmlly obtained in practice . 

FigTire 2 also shoves the advantage of using thin airfoil sections 
to absorb higher power. The power loading at =0.8 is 5l8 horse- 
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power per square foot for a 6-percent— thick airfoil section and 
390 horsepower per square foot for a 15-percent-thlck airfoil 
section. This conclusion regarding the effect of thicknoss applies 
only for airfoils operating at the design lift coefficient. 

Because of high peak pressures at other than the design lift coeffi- 
cient, the thin airfoil sections may not iiave the advantage shown 
hy figui'e 2. 

Conventional airfoil sections .— As pointed oi’t hy Heaslet the 

idealized sections are not practical airfoils, chiefly hocause the 
elliptical thickness dlstrihution would cause flow separation on 
the afterbody; thus, high drag would result. Tiie idealized section 
of a given thicloiess ratio, furthermore, repi’esents a different shape 
airfoil for each lift coefficient. Thus, the curves in figure 2 
for 6-percent thickness represent not one airfoil shape htit an 
infinite number of airfoil shapes of 6-percent thickness. In order 
to determine how near conventional airfoils approach the ideal! ze<i 
section in critical Mach nxraiber and power loading, tliree conventional 
airfoils are compared with the idealized sections of the same thickness. 

The three conventional airfoils are a Clark Y airfoil section 
having 11.7— percent thickness and 3*5— percent camher, a Clark YM 
airfoil section with the same thickness dlstrihution hut with 
9— percent camher, and an NACA l£^serlos airfoil of 12— percent 
thickness, 5*5— percent camher, and a design lift of 1.0. Tlie airfoil 
shapes are shown in figure 3> The coordinates of the first end second 
airfoils are given in reference 3 end the coordinates of the last 
airfoil are in reference k. 

The critical Kaoh number as detonnined theoretically is plotted 
in figujre k for the conventional airfoils and the 12— percent-thick 
idealized airfoil. The critical Mach number for tlie idealized— airfoil 
curve decreases with an increase in lift coefficient. Kach conventional 
airfoil is seen to approach the cui-ve of the idealized airfoil over 
a certain range of This favorable range coincides approximately 

with the design or ideal lift of the peirtlcular airfoil. At low 
values of the lift coefficient the critical Mach numbers may he low 
because of t-’ie high velocity at t’le front lower surface of the 
airfoil. 

Figure 5 gives the maximum power loading for the same airfoils. 

The power loadings for the high— camher airfoils become favorable 
only at the higher lift coefficients. At the high lift coefficients 
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the hlgh'cainber sectionB hava Tiller power loadiagB .toaii the low- 
camber section. Some doubt exists, however, as to the maociraura lift 
coefficient that can be used in practice without separation of flow 
resulting in low efficiency. 

The design of the whole blade for a given operating condition is 
beyond the scope of the present paper. If the entire area of a 
propeller or fan blade is considered, the average power loading is 
much lower than the maximum power loading which refers only to the 
loading of an Isolated blade section operating at its critical 
speed. This difference occurs because the inboard sections are 
operating at much lower stream Mach numbers and because the helix 
angle varies along the blade. The average power loading is probably 
increased by reducing the tip lift coefficient and increasing the 
tip Mach number j thereby, the power loading of the inboard sections 
is increased. 


Maximum Power Loading as Limited by Flutter Speed 

Flutter is a self -excited oscillation of a body caused by enorQr 
absorbed from the air sti'ea:'!!. This oscillation is usually very violent 
nnri destructive . The two principal types of flutter are classical 
flutter and stall flutter. Classical flutter is an oscillatory 
instability of an airfoil operating in a potential flow. In general, 
such flutter requires at least two coupled degrees of freedom, such 
as bending and torsion. Stall flutter is caused by separation of 
flow and occurs on airfoils operating near or in the stall condition 
of flow. This type of flutter requires only one degree of freedom, 
usually torsion, and is generally attributed to the hysteresis in 
the lift curve near stall. 

In reference 1 the classical flutter speed and the divergence 
speed are shown to be almost the same for propellers and fans. A 
design sufficiently rigid to preclude divergence is usually safe 
against classical flutter. The aerodynamic moment may however 
appreciably twist the blade at much lower speeds. This twisting may 
change the angle of attack sufficiently to cause stall flutter. 

In reference 1 the stall flutter was seen to occur on a conventional 
Clark Y airfoil section whan the blade twisted sufficiently to 
increase the lift coefficient to about = 1-0. This value cannot 

be talcen as an absolute limit since the stall characteristic depends 
on RejTcolds nutaber, Mach number, and type of airfoil. No aerodynamic 
twisting raoiaent would be present on the blade if the aerodynaxiilc 
center of pressure of the blade section coincided with the center of 
gravity of the section (reference 1) . For this condition the blade 
does not twist until the divergence speed is reached. The relation 
for this lift coefficient of no twist is given in reference 1 as 
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(B) 


vhers 

Cb’^/ 4 pitching-moment coefficient about quarter -chcrd point 

Xjj.g. location of center of gravity as measured from leading edge , 
chox'ds 


If the lift coefficient is greatei- than C; 




the aii’foii. flutters 


at redx’ced speed with positive stall; if the lift coefficient is less 
than flutter occurs at reduced speed with negati've stall. 

In either case, the maximum power that the blade absorbs occurs when 
the blade is operating at its lift coefficient of.no twist C]^ . 

The following equations, based on equations from refei’ence 1, 

have been used to calculate the lift coefficient at flutter speed 

for any value of Cm : 

u 



- C' 



(9) 


and 


= Cl 


^iv 


(Cl - Clu^^ 


(10) 


where 

lift coefficient 

untwisted or design value of lift coefficient 

Ct lift coefficient for ideal no-twist condition 

q/q^i^ dynamic pressure of operating speed divided by dynamic 
pressui'^ of divergonce speed 
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From equation (1) the power loading is seen to he proportional 
to Since q is proportional to the power loadin|^or 

a given value of q,.„ is therefore proportional to 

Figure 6(a) gives the relative maximum power loading at flutter speed 
for the condition that is greater than Tlie curves 

are calculated hy use of equation (lO) on the assumption that flutter 
occurs when Cl = 1.0. As -i/aa pointed out previously, all airfoils 

do not flutter at exactly Cl = 1.0. The use in this analysis of 
this value for Cl as the lift coefficient at which flutter occurs 
is for purposes of estimating the vai’lation of pov/er loading for 


various values of 


CLu ““i =Lu, 


is less than Ct 


the operating lift coefficient 


decreases with increased speed and the maximum power loading drops 
off very rapidly as Cp^ differs from Cq^^. In this case the 

maximum power does not occur at the flutter speed hut at some point 
helow the flutter speed. The maximum power loading for various 
values of Ct is plotted in figure 6(h). The curves in this 

■^Ut - 

figiare were obtained hy graphical methods and use of equation (9) • 

Figure 6(a) shows that if a hlade section having a value of 
Cq^ = 0.4 is operated at Cq^ = 0.8 the hlade ahsorhs only 0.2 of 

the power it would ahsorh if it were operated at Cp^ = 0.4. Figure 6(h) 
shows that If the hlade section having a value of Cp^^ =0.4 wore 

operated at Cp^ = 0.2 the hlade ahsorhs only 0.02 of the power it 
ahsorhs when operated e.t Cq^ = 0.4. This example Illustrates that 
for a hlade to ahsorh maxiiiumi power, Cp should equal Cp^ j and, 
if a deviation exists, the value of Cp^ should he greater than Cq^ 

rather than smaller. These conditions are in line with standard 
practice . 


The dynamic pressure of divergence speed is that pressure 

at which the aerodynamic -moiTient stiffness of the hlade equals the 
torsional stlfx''ness at the hlade and may he obtained from vibration 
date as indicated in references 1 and 5* Neforence 1 gives an 
expression for the divergence speed in terms of the torsional 
stiffness of the hlade as follows; 
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^dlv " 


K 


Lc2 


<i£L 

da 


(xc, 


e- 


I') 

hj 


( 11 ) 


where 

L representative length of hlade^ feet 

c chord, feet 

^Q.Q. location of center of gravity as measured from leading 
edge , chords 

K torsional stiffness of "blade, foot-pound per I’adian 

The torsional stiffness K is proportional to where G is 

the shear modulus of elasticity of the material and J is the torsion 
modulus of the section. For thin sections, J is approximately 

O 

proportional to t-^c where t is the thlclaiess of the section. 
Equation (ll) may "be written 



( 12 ) 


Equation (12) shows that varies directly as the shear 

modulus of the mterial., as the cu"be of -the thickness ratio, and as 
the square of the chord-length ratio and varies Inversely as the 
distance of the section center of gravity from the quarter-chord 
position. For example, if the hlade thickness is increased from 

6 percent to 15 percent, qd.iv is increased "by a factor of 

or 15-6. If a "blade is made with a sufficiently large thlclaiess 
ratio or sufficiently large chord -length ratio, the flutter problem 
can be eliminated entirely; This inci^ase of ratios, however, is 
done at the expense of more weight or of poor aorod^mamic character- 
istics and lower critical speeds. 

Tlie stall-flutter speed expressed in terms of dynamic presaiore 
may be obtained from equation (lO) on the assumption that flutter 
occijTS at a given lift coefficient The equation is 
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For positive stall flutter, If Cj^ 

equals C, the term in the parentheses equals unity, which is 

•*-'Uj 

its maximum value. Any deviation from this condition reduced the 
flutter speed. By proper design of the chordwiso load distrihution, 
may he made identical to the ideal or desigi lift coefficient 

of the hlade section. The ideal lift corresponds to the lift at the 
ideal angle of attach, which is defined as the angle at which the 
front stagnation point is at tiie leading edge of the airfoil. If 
the section is operated at this ideal lift, the liigjiest critical 
speed and "best flow conditions are obtained for that operating 
condition. (Sea reference 6.) 


For the homogeneous conventional airfoil sections considered 
in the present paper, Cjj^ is almost equal to Cj • Table I 

- i -Uj 

gives the values of Ct o,nd Ct for the three airfoil sections pre— 

■“Uj ■‘■'I 

viously discussed. The values of C-^ are calculated by equation (o). 

The values of C«, and Ct were obtained from potential 
“c/4- 

calculations. As seen from table I, operating the airfoils at 
has no disadvantages since Ct and Ct are a.lmost the sane. 

-Ail 


The value of q^^^ that should be used depends on the 
application. In general, a very higli value of q^^^ gives a 

etlffer blade having a longer operating life, particularly in 
applications wliere vibration is severe. A higli value of qj^^y 

la also necessary whore the operating conditions varj' widely, that 
is. where tlie value of Ct varies over wide limits. Equation (13) 

shows tliat for such conditions qp may be only a small free ti on 

of q-,. . In some wind-tunnel fans end axial -flow compressors, 

Mlv 

ho'^TOver, the operating condition with respect to the blade lift 
coefficients does not vary appreciably. In such cases a lighter 
blade having a lower value of q^^y iJia-y be successfully used. It 

appears that the most economical and efficient design would be one 
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in which 



of the "blade wei-e cquEtl to the rainimuin operating 


lift coefficient for maximuiii speed. In general, such a desigri gives 
the highest critical speed as well as the greatest flutter safety 
margin foi- a given value of 


SXPERDiENTAL STUDIES OF FAN FLUTTER AT EECxH LOADING 
Equipment and Test Procedure 


Wind-tunnel-fan Models were tested in tho same open-return 
tunnel used for the tests of reference 1. A diagi*aaiatic shetch of 
the test setup is shown in figure ?• The "blade tips were illuminated 
"by stroboscopic light and observed in operation b^r means of a 
window in the tunnel \rall. Tho blade load and Ct of the blade 

were varied by me.ans of slats in the tvumel exit. These slats 
change the tunnel velocity and thei*eby change the angle of attack 
of the fan blades. 


The fans are i^iade of laminated Sitka spruce and have a specific 
gravity of about O .5 and a diaiaetor of 45 inches. B3.adG-forM curves 
for the fans tested are shown in figure 8. Tho fans are as follows: 

Fan A is a six -blade fan having conventional Clark Y airfoil 
sections. This faia is tlie same as tlie fan dealgna-ted in reference 1 
as propeller A. 


Fan B is a four-bleB.e fen having Clark "JfM airfoil sections of 
9-p©rcent camber. This fan has the sane ble.de dimensions and lift 
distribution along the blade as fan A. Fan B is also operated as 
a six -blade fan to obtain seme flutter points at low lift coefficients. 
All data reported are talcen for the four -blade fan except one datura 


point at 



in figiue 0 . 


Fan C is the saxae four-blade higli -camber fan as fan B but with 
the difference that one blade was weakened by cutting out part of the 
spruce and Inlajdng a balsa insert. The cut-out was necessary to 
weaJeen the blade in order to obtain fluttei* over a vide reuige of 

The flutter speeds of fan C refer only to the flutter of the weakened 
blade. Cutting and inlaying the blad,e warped the blade sufficiently 
to give it- about 11 -percent camber. The value of for this 

section is therefore not the scima as that given in table I for the 
9-percent-cairber Clark Y14 airfoil section. 
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Results 

Comparativ© j^ieasurenents were aiade with the six -"blade fan A 
and. th.6 four -Had© fan B to determine tiie flatter cliaracteristicSj 
maximum power loading, maximimi lift coefficient, and efficiency of 
the two fans. Tlie power loading is expressed in teimia of horse- 
power per square foot. Each "blade has an area of about 1/4 square 
foot. 


The results of these tests are given in figures 10(a) to 10(e). 

The efficiency is defined as the useful power divided by the motor 
output as measured by a strain-ga{^ dynaiaon^ter . The^power is the 
Integral of dQ AP.? where Q is the quantity rate 01 air flow througli 
the fan and Ap is the pressure rise through the fan. The absolute 
values of these measui’ements are not very good, but the same 
systematic errors apply to both fansj consequently, the ef.-. iciencies 
shoHd be useful only for comparison purposes. The f^ losses are 
due to the drag of the airfoil sections and tne rotatxona.1 energy' 
loss. Since the airfoil sections were small and rough (maximum 
Reynolds number of 1,000,000), the section drag losses were proba"bly 
hi^er than for full-scale fans. From the consideration of efficiency, 
there seems to be little choice between the two fans. 

The lift coefficients are given for the 0.9 radius. The lift 
coefficients are obtained from total -pres sure measurements tahen 
behind the fans and by use of equation (6) . The rotational- 
■yelocity interference factor a' is low for the tip section of 
the fans tested. This factor was therefore neglected in the 
determination of W for the calculation of the lift coefficient 
(a' =0). 

Since the model Reynolds nuvd)era were small and the models were 
rough, the maxiraun lift coefficients are probably less than wouild be 
obtained with larger Reynolds numbers and s^iToother models. F^. A, 
having a conventional Clarlc Y airfoil section, gave maximum Hit 
coefficients of 0.8 at the 0.9 radius and of 0.9 on some of the 
inboard sections. Fan B, having the 9“psrcent-cariibor Clark "Yin 
airfoil sections, gave maxiimim lift coefficients of about 1.0 ab 
the 0.9 radius and of 1.4 at some of the inboard sections. Although 
the higher -camber fan does give somewhat higher maximum lifts, 
these high lifts are obtained at some losses in efficiency. The 
prediction of how much lift would be carried by fans operating at 
high Reynolds numbers, without separation of xlow occurring wxoh 
a resulting loss of efficiency, is dlfficulLt b'ccause of some uncertainty 
of the effect of Reynolds number on the maximum lift. 

The blade -tip twisting as a function of the lift coefficient 
for constant fan speed is given for fan E in figure 11. The lift 
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coefficient is £;iven for the O .9 radius and is calculated from total 
pressure EisaourGinents as pi'eviously discussed. The twist at the 
tip was ohserved hy means of a strohoocope and was moasured w?th 
a protractor on a telescope. Figure 11 shows that for constant 
operating speed the Hade twisting deforiiiation Increases the ivagl® 
of attack at the high lift coefficients and decreases the angle of 
attack at the low lift coefficients. Note that a lift coefficient 
exists for which the Hade twist is zero. This value of is ly 

definition the experimental lift coefficient of no twist Cj_^ . 

The value of Cj^ for fan B is approxtuately 0.8 (fig. ll). This 

value is less than the theoretical value, I.I 6 , for the Clark 
airfoil section given in taHe I. This discrepancy is "believed to 
"be due to section boiandary -layer effects- The experimental vHue of 
for fan A is given in reference 1 as 0.37- Although f an B 

has an experiirental value of Ct leas thtn that predicted frora 

theorj>-^ the experimental value of Cy is about twice that for 

fan A. A f;an with the same section as fan B can be erirployed at high 
lift coefficients with a greater margin of safety than can a.low' 
caiabei’ fan section. 

The decrease of with tip Mach number as shown in figures 10 ( 

and 10(b) for fan B is obviously caused by the blade fcristing and 
decreasing the angle of attack because the fan is operated^ at a lift 
coefficient below the eroorimentally determined value Ct =0.8. 

The decrease of at higii Mach numbers shown for fan B in 

figure! . 10(c) is primarily a compressibility effect augnentod by the 
blade twisting in a negative direction because the center of pressure 
has moved back at the high Mach numbers. Figures 10(d) and 10(e) 
give resHts for both fans partly stalled witja resulting low 
efficiency. 


The power loading of the fans as given in figure 10 Is expressed 
in terms of horsepower per sq^uare foot. The maximum values obtained 
in thesS te.sts are much below those shewn in figures S and 5 which 
give theoretical values calculated on a basis of a tf-5^ helix fuigle 
and the critical Mach nuiaber of an isolated section. For a fsn of 
which the Mach number and helix angle vary greatly with the radius, 
the maximum values given in figures 2 and 5 c-onnot be attained. For 
an axial-flow compressor, however, of which the blades are short 
and made to operate' at' easentiall;^ a cous'taht helix angle and Mach 
number-, it my be possible to attain a power loading approaching 
that given in figures 2 and 5- . ■ 
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The critical Mach numbers obtained from figure 4 are gi.ven in 
figure 10. The values of used for obtaining are taken 

at Mj. = 0.4. Since changed with tunnel conditions, the 

critical Mach numbers are different. The critical Mach numbers of 
both the upper and lover surfaces of fan B are also Indicated. For 
fan B at low lift coefficients a high velocity peak occtxrs on the 
lower leading-edge surface] thus, a low critical Mach number results. 
This localized supersonic region does not seem to have much 
significance. The conclusion is in agreement with the conclusion 
in reference 4, which states that the Mach numbers at which large 
clianges in airfoil cliaracteri sties occur are difficult to predict 
especially when sliarp pressure peaks exist at the leading edge. In 
figure 10(b) no significant changes in the lift coefficient or the 
efficiency are seen to occur until the critical speed on the upper 
surface has been exceeded. This drop in efficiency is probabl5’’ due 
to flow separation caused by shock, as po5.ntod out in reference 7* 
These tests indicate that the critical speed of the upper surface 
as calculated by two— d-imenslonal theory is essentially a limiting 
speed for efficient operation of ducted fans. 

One of the most significant results of the investigation is 
given in fig\ire 9. This figure gives the maximum power absorbed in 
terms of horsepower per square foot for fans A and B as a function 
of the lift coefficient at 0.9 radius. The curve for fan A 

represents the power absorbed as limited by flutter . The maximum 
power loading is about 100 horsepower per square foot of blade area 
and occurs at a value of of about 0.47- This result is in 

fair agreement with the theory which shows that the maximum power 
loading should occur at about 0.37 for fan A. Fan B did 

not flutter at values of between 0.5 to 0.85] consequently, 

this part of the curve is shown as a short-dash line. In this 
range the power is limited by the maximum speed of the motor. This 
maximum power loading for fan B is beyond the limits of efficient 
operation because the fan is operating in the supercritical speed 
region. The flutter points were obtained at values of Cp^ of 0.4, 

0.88, and 0.92 and this part of the curve is shown as a solid line. 
The long— dash curve gives the power loading at the point where the 
fan efficiency has dropped to 80 percent. The power loading of 
fan B reaches a maximum value of 150 horsepower per square foot 
which is considerably greater than that reached for fan A. 

Tlie flutter point for fan B at Cp^ = 0.4 was obtained with 
a six— blade fan] the points at Cp^ = 0.88 and = 0.91 were 

obtained with the four-blade fan B. The steep slope of the power- 
loading curve for fan B at the high and low lift coefficients 
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indicates that a fan having hish-caniber sections can he used 
successfully to absorb a large amount of power over a certain range 
of hi^i lift coefficients but that such blades have undesirable 
flutter characteristics at the low lift coefficients. 

The vibration frequencies and flutter speeds are given in 
table II. The flutter speeds are calculated by the method of 
reference 1 when the 0.8-radius station is used as the reference 
station. The maximum ani minimum measui'ed flutter speeds are also 
tabiilated. The mlninm speeds are those obtained with the blsides 
completely stalled. 

The maximum flutter speed for fan A is somewhat less than the 
calculated classical -flutter speed corrected for compressibility. 

Fan C checks the calculations quite closely. Fan B had some 
Incipient flutter at a tip Mach nmber of O.67 when operating with 
the tunnel open. This Mach n'umbor corresponds to tfie approximate 
critical Mach number of the upper blade surface for tunnel-open 
condition. As the speed was increased into the .supercritical 
region, this flutter disappeared and the blades operated very 
smoothly to the top speed of the motor which cori’esponds to a tip 
Mach number of 0.86. At these top speeds the blade lost most of 
the lift near the tip and, also, tv/’isted in a negative direction. 
Although the ducted fan operated smoothly in the supercritical 
region at speeds above the classical -flutter speed, it does not 
necessarily follow that such will be the case for tlifee -dimensional 
bodies. 

The last row of table II gives the minimum measured stall -flutter 
tip Mach niimber divided by the calculated classical-flutter tip 
Mf . - 

Mach numbei’ — . This number is of particular interest when 

Mf 

^cal 

fans are considered which Eiay be required to operate in the 
completely stalled condition. No reliable theory is available, 
unfortunately, for predicting this minimum flutter speed, which 
veries for different blade -section shapes, blade plan form, blade 
materials, and so forth. At present any design which is intended 
to operate in the stalled condition should be tested by overspeed 
whirl tests . 


CONCLUSIONS 


The results of the analysis and of tests made of t^/o \find-tunnel- 
fan models to determine the efi’ect of flutter speed and critical 
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Mach lUBnhor on ths poi7©r loading of the hlade sections indicated 
the following results: 

1. In spite of the fact that hi^ier lift coefficients result in 
lower critical speeds, the power that a hlade section ahsorhs at 
the critical speed is a riaxinum at the liig’iest lift coefficient the 
hlade section can develop. This fact apjjlies to idealized sections 
and to conventional sections operating near the ideal li.ft. 

2. If flutter limits the operating speed, the fan ahsorbs the 
maximum power when the "blade is designed to operate at the lift 
coefficient of no twist. When this lift coefficient and tlie 
ideal or desj.fii lift coefficient are identical, a "blade operating 
close to its ideal lift has a riaximuia critical speed as well as a 
maximum flutter speed. Any material deviation from the ideal lift 
coefficient resul.ts in a greatly reduced maxiinitm power loading. 

3- The anticipated increase of power leading of the high- 
camber "blade was o"btai.ned in tie tests. There was little difference 
In the efficiency between the ClaLrk Y and the hlgli-camber Clark "IfM 
fans, but the hl£^i-ca:vLber blades developed somewhat higJ'ier ira:oimura 
lift coefficients. The lilsh-camber blades, how'ever, were found to 
have very poor flutter characteris,tics at the low lift coefficients. 

4. The fan efficiency decreased rapidly after the sonic velocity 
on the upper blade surface of the fan sections was exceeded. 


Langley Memorial Aeronautical Laboratory- 

National A*^"vlsory Committee for Aeronautics 
Langley Field, Ya. , April 16, 1947 
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TABLE I . - COEFFICIENTS AID LOCATION OF CENTER OF GRAVITY 
FOR AIRFOIL SECTIONS 


Airfoil 

Cnic /4 

^c.g. 

0 

}-I 


Standard Clark Y 

-0.085 

0.44 

0.50 

0.45 

Clark YM 

-.22 

.44 

1.29 

1.16 

NACA 16-1012 
_ 1 

-.25 

.43 

1.00 

1.06 


TABLE II.- VIBRATION FREQUENCIES AJIB FLUTT'ER SPEERS 


Items 

Fan A 

Fan B 

Fan C 

First tending froquoncy, ops 


6)+ 

76 

Second tending freq^uency, cpa 

246 

245 

220 

First torsion frequency, cps 
Calculated classical -flutter or 

355 

365 

287 

divergence speed at 0.6-radius 
station, fps 

7T2 

795 

624 

Calculated classical -flutter tip 




Mach numter, Mf , 

' cal 

0.86 

0.89 

0.70 

Calculatod classical -flutter tip 




Mach numter, corrected for 
compreasitility 

0.76 

0.73 

0.64 

MaxiBium measiired flutter tip 




Mach nvcntei' , Mp 

’ -tjnax 

0.71 

<a) 

0.66 

Minimun measured stall -flutter 




tip Mach numter, M;- 

0.34 

0.44 

C.28 

’ -"min 



^^min 

^f 1 
cal 

o. 4 o 

0.50 

o. 4 o 


^■No flutter up to maxlnum operating speed of tip Mach numher 0.85. 
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Figure 1.- Diagram of velocities and forces considered In 

analysis. 
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Figure 2.- Maxlmvun theoretical power loading and 
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Fig. 3 
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Figure 5*“ Shape of airfoil sections 
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Fig. 5 
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Figure 6.- Concluded. 
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Flgvire 7 «“ Cross-sectional view of wind tunnel used for 
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Fig. 8 



0 ) 


CCL 


Figure 8.- Blade-form ctirves 
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Tunnel blocking, 35 percent. For fan A, M(.p - 0.68; for fan B 
(upper surface), = 0.6ij.; and for fan C (lower surface), 

Her = 


Figure 10,- Continued 


NACA TN No. 1330 


Fig. 10c 
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(c) Timnel blocking, percent. For fan A, M^p = 0,66j for fan B 
(upper surface), M^p = 0.6l; and for fan C (lower surface), 
Mqp “ 0 • 50 . 


Figure 10,- Continued 
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(upper surface), = 0.62; and for fan C (lower surface), 

Mcr = 0*^9. 


Figure 10,- Continued 
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